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j SUMMARY 

A numerical calculation has been carried out to evaluate the 3x3 cross-section 
matrix involved in the electron impact excitation of the ground state of H atom to 
the 2s and 2p levels. The method of solution is that of atomic eigenstates expan- 
sion. In this paper, instead of the iterative technique used by other authors, 
the definite integral terms in the coupled radial differential equations are 
eliminated through some linear transformation of the radial functions, thus avoiding 
iteration of these equations. The accuracy of ,the numerical integration is tested by 
satisfying the equation of reciprocity and the equation of continuity of currents with 
an error-to-value ratio less than 1 per 1000 on the average; and the maximum of 
this ratio, except for a few cases, has been kept below 5 percent. The results are 
in agreement with the results of an iterative technique. 

To evaluate the effect of the long-range and the centrifugal potential, a simple 
perturbation theory is developed. The six cross sections ls-2s,ls-2p, is -is, 
2s - 2s , 2s -2p, and2p-2p are tabulated. The2p-2p cross section requires the 
solution of the sets of differential equations with different parities. With the validity 
of the eigenstates expansion assumed, it is found by comparison with the eigenstates 
expansion calculation that the Born approximation, despite its simplicity, gives 
meaningful results for low and close-to-the-threshold energies of the bombarding 
electrons. The effect of the exchange potentials on the cross sections is also in- 
vestigated. Finally, an interesting structure of the Is -2s excitation cross section 
above threshold is found. 
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THE HYDROGEN ATOM INDUCED BY ELECTRON IMPACT* 


by 

Kazem Omidvar 
Goddard Space Flight Center 


INTRODUCTION 

Calculation of the excitation cross sections in atomic hydrogen by electron impact corre- 
sponds to the solution of the problem of three interacting bodies: one proton, and two electrons. 
By taking the position of the proton as the center of mass, the problem will reduce to the task of 
finding the nonseparable wave function of the system of the two electrons with an attractive cen- 
ter of force. Such a solution has not been found. However, if this wave function is expanded in 
terms of the eigenstates of the hydrogen atom, the coefficients of the expansion, which are func- 
tions of the position vector of the free electron, can be found through numerical integration. When 
an infinite number of terms is included in the expansion, the solution to the problem is exact. 
Furthermore, the expansion has the advantage that the asymptotic form of its coefficients is 
automatically the asymptotic form of the free electron wave function scattered from different 
atomic states, which are simply related to the excitation cross sections. 

In this paper, atomic states is, 2s ,2p are included in the expansion and, by antisymmetrizing 
the two electron wave functions according to the exclusion principle, some contribution from the 
continuum in the expansion is also taken into account. The first calculation of this type was per- 
formed by Marriot (Reference 1), whose expansion consisted of the is and the 2s states in order 
to calculate the is -2s transition cross section. This calculation was extended by Smith (Refer- 
ence 2) to higher total orbital angular momenta of the system. Percival and Seaton (Reference 3) 
have formulated the eigenstate expansion technique in general and have tabulated the coefficients 
of the integro- differential equations for s , p, and d atomic electrons. Burke, Smith, and Schey 
(References 4 and 5 f ), using the equations of Percival and Seaton for the three states is, 2s , 2p, 
have integrated the resulting integro- differential equations. In this paper we solve the same dif- 
ferential equations by a linear transformation of the differential equations in order to avoid the 
need for iteration of these equations (Reference 6). $ 


•Also has appeared in a condensed form in the Physical Review , Vol. 133, Feb. 17, 1964. 

tA similar calculation has been performed by R. Damburg and R. Peterkop; this will appear in the USSR Journal of Experimental and 
Theoretical Physics. A different method to calculate the ls-2s electron impact transition cross section in hydrogen is being con- 
sidered by L. Kyle and A. Temkin, adopting the nonadiabatic theory of electron scattering developed by A. Temkin (see References 
21, 22); the calculation is in progress. 

tA similar calculation has been carried out in Reference 6. Here the L = 0,1 cases have been solved by noniterative, and all other 
cases by iterative, methods. 
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The numerical integrations were carried out for all partial waves, while in higher partial 
waves the Born approximation was used. The transition between the eigenstates expansion cal- 
culation and the Born approximation takes place when the results of the two calculations agree 
closely. 


FORMULATION 

Derivation of the Differential Equations 

Since spin orbit interaction of the electrons is neglected, the total orbital angular momentum 
L and the total spin angular momentum s are separately conserved. We can then divide the inter- 
actions into antiparallel spin states, where S = 0, and parallel spin states, where S = 1 . We deal 
with spatial wave functions of the electrons only, and for brevity we call the orbital angular 
momentum the angular momentum . 

Neglecting the motion of the proton of the hydrogen atom and taking its position as the origin 
of the coordinate system, the Schroedinger equation for the system can be written 


[H-E] 0( r lf r 2 ) - 0 , (1) 

where and r 2 are the position vectors of the bound and free electrons; and in atomic units 


H-E 






( 2 ) 


where E is the total energy of the system and r I2 is the distance between the two electrons. We 
expand the total wave function \p(r t , r 2 ) in terms of the eigenfunctions of the total angular 
momentum L , 


H* i- r 2 ) 


l>L ( r l- r 2 ) • 

L=0 


(3) 


Since these eigenfunctions are orthogonal and distinct, substitution of Equation 3 in Equation 1 
gives 


[H-E] 0 L (r lf r 2 ) = 0. 


(4) 
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The explicit form of (r,, r 2 ) is given by 


( r i* r a) - 


n l<^ 2 


c^4>A 


r l) 


m l m 2 


X r 2 ->u( knj ^,r 2 )F V2 (n 2 ) , 
fj) = r' 1 r,) (n,) . 


(5) 

( 6 ) 


Here ^(rq mj , rj) is the hydrogen atom wave function with radial part rj' 1 P[n l <£,, t 2 ) and angu- 
lar part Y (n 2 ) and quantum numbers n-l 1 mi ; r " 1 u^k^ , r 2 j is the radial part and m (n 2 ) is 
the angular part of the free electron wave function with quantum numbers k n l 2 m 2 . The relation 
between the wave number k n and n, is given by 


k 



(7) 


^ L 

Finally the constants * = ( / ^ l ^ 2 m i m 2 1^** )> w*** 1 M representing the total magnetic quantum number, 
are vector coupling coefficients which make the linear combination of the products of the one elec- 
tron wave function in Equation 5 the eigenfunction of L. In the problem under consideration, 
nj = 1, 2 ; / i l - 0, l; 4 2 = |l- "fcj | , , + ; m x u - <i v and m 2 = - * * * > ^ 2 * 

To make the total wave function symmetric for antiparallel spins or antisymmetric for parallel 
spins, the operator P 12 interchanges r 1 and r 2 while /3 is +1 for the first case and is -1 for the 
second. 


By taking L perpendicular to the z-axis, M = 0 andm 2 = -m x . Equation 5 can then be written 


( F 1 ' r 2 ) (* + ^12) 



c vv*( n i <t i m i- r 0 




* r 2 _lu ( k n I /e 2- T i) Y V 2 fa) * (8) 

In order that ^(r, , r 2 ) closely approximates the exact wave function, we minimize the ex- 
pectation value of the energy operator with respect to the radial parts of the free electron wave 
functions, 


5 J*^l ( r i- r z) [H_E] '/' L ( r i> r 2 ) d3p i d3 


= 0 


(9) 
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It has been shown by Kohn (Reference 7) that the differences between the scattering amplitudes 
obtained from these equations and the exact scattering amplitudes are quadratic in the difference 
between^ (r 1 , r 2 ) and the exact wave function. When the variation is carried out inside the in- 
tegral, we obtain 



">!• r i) *4 a ( n 2 ) [H " E], M r i* *- 2 )d3 Pl dn 2 


0 . 


( 10 ) 


By means of Equations 2 and 8, the Schroedinger equation for the hydrogen atom, 




-Cj 


r l) " — ; m,, r,) 


and Equation 7, Equation 10 reduces to 


(id 



Q ^1^2 L 
mj-nijO 






'i)v (n 2 ) ( i + /3p 12 ) 




x 



V (V + 1) 




F l) 


x r a -> u(k n; i ' , r 2 ) Y^. (fi 2 ) r, dfi 2 = 0 , (12) 

where V 2 is the radial part of V 2 By orthogonality of the hydrogen atom and spherical harmonics 

T 2 

wave functions, the relation (Reference 8) 



the integration by parts of the exchange terms, and the relation* 


C ^V 

“irijiTjO 


l-v* 2 

v ) m.-m. 0 


(13) 


( 14 ) 


’Reference 8, Equation 3.16b. 
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Equation 12 leads to 


^2 (^2 + 2 
V 2 - + k 2 + — 


** (n,-^ m,, r,) y/ (fi 2 ) 


H U (V” r a ) 


- 2 



c "V2 l c W l 

nij-mjO 


vW »i' 


+ & 


U ( k ^J' r 2 ) 

r,) ~ f" Y^(n 2 ) 

u ( k n' "V- r l) 

+ ^( n l' •£/ m,' , P 2 ) X J- \ i , m . (n,) 


2^(-) L ^ 2S (' e +k nj) 

■vW 


d 3 Tj dH 2 


J o p ( n i*i' r i) r 2 ) x u( k „; V. r x ) dr, = 0. (15) 

If 1/rjj is expanded in terms of the Legendre polynomials and use is made of the addition theorem,* 
we obtain 


a> 



\=0 



M n .) C (n 2 ) 


(16) 


In this expression e n is the angle between the position vectors r, and r 2 at the origin, and r< is 
the smaller and r> is the larger of |r, J and jr 2 | . We also introduce 


f r 2 

(iv£ n' •£ ' , r 2 ) = r 2 ” ( ^ +1 i I P(rrf, P(n' ■V , Tj) dr t 


+ r 2 ^ j* P(rvt, r ^ P(n' V , r x dr x . 
J *2 


(17) 


^Reference 8, Equation 4.28. 
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Then it follows that 


f” P ( ni > r i) p ( n ' ' t ‘* r,) 

00 +K 

L £ y v («. ) 

J r i 2 d 1 


X= 0 m- ”X. 

By means of Equation 18, the relation* 

r * 

"(2^ + 1) (2-tj + l)~ 

J \. a Y V 2 Y V, dQ = 

4tt (2-e 3 + i) 


(18) 


1/2 


c Ws c Ws 

m l m 2 m 3 OOO 


(19) 


and the definition 


/•CD 

^iv£|k n = J P(rv£, r ) u^k n / fc # , dr , 


( 20 ) 


Equation 15 when multiplied by r 2 gives 


d 2 l l ( l 2 + l ) 


dr 2 


+ k 2 + 


J ( k n,^2- r 2 ) 


- 2 


2-t 2 + lY 2 


2l 1 + 1 


00 +\ 



U 2 L tjtj L 

^mj _ m|0 m i _m /o 


_l L I 

m \ m t ^ = 0 fx~-X. 


^2^+1 V ^ ^ XV; u 2 t' , . W|f . , \ 

i ' + j I V.;., C 000 %», C 000 \ K n l *1 - r 2 ) • r 2 ) 


2-t,' + l\ 1/2 


+ b\ — - ) r^ 2 ^ 1 c^ 2 ^ 1 c^i' pi ‘ x • \ lx v 

2-tj' + 1/ ^"V”! C »00 C Wm 2 m,' C 000 P ( n l ^1 • r 2 ) \ [ n l ^1 k nj' -tj , r 2 


+ P (-) L ^ 8(4,' -t 2 ' , -t 2 -fj) ^ +k n 2 ,j />(nj' -t/ , r 2 ) (n, l x |k^ 1 2 ') - 0 

"l'W 


( 21 ) 


In the exchange integrals above we have defined /^k nj -f 2 , r j as 1 2 , 


•Reference 8, Equation 4.34. 
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The vectorial equations 


L = l x + l t 


( 22 ) 


where L is constant but l x and l 2 take the values given before, can be divided into two groups, one 
with L - l x - l 2 even and the other with L -if l 2 odd. Since the total spatial wave function has the 
parity of + 1 2 , in the first group the wave function has the parity of L and in the second a parity 
opposite to L . By conservation of parity we have two distinct groups of interactions. In this 
problem, where is, 2s, and 2p states of atomic hydrogen are taken into account, it is easy to see 
that, when L - i t - l 2 is even, the set of quantum numbers l 2 has four values: one for each of 
the Is and 2s states, and two for the 2p state. WhenL-^ -l 2 is odd, k n l 2 has one value which 
corresponds to the elastic scattering of electrons by the 2p state of the hydrogen atom. Equation 
21 is evaluated for these cases, and the resulting differential equations are listed in Appendix A. 

In evaluating Equation 21, it should be noted that the c m 2 m 3 coefficients are subject to the condi- 
tion that l 2 l 3 form a closed triangle and m 3 = m 1 +m 2 .* This limits the summation over k and 
considerably to few terms only. Summation over m 1 , m/ , and & is carried out using the numer- 
ical values of the C coefficients given by Condon and Shortley (Reference 9). 

Percival and Seaton (Reference 3) have derived the same differential equations for the scat- 
tering of free electrons by atomic s , p, and d electrons in the hydrogen atom using the theory of 
irreducible tensor operators to evaluate the interaction terms between the two electrons in the 
differential equations. The calculation becomes considerably simpler in this way. The results 
of the two methods are identical. 


In the rest of the paper, except the section on page 14, we discuss the solutions to the four 
coupled differential equations given in Appendix A and which arise when L - - 1 2 is even. The 

single differential equations for odd are derived in the excepted section (page 14). Its 

numerical integration can be treated as a special case of the four coupled differential equations. 

When the integrals representing the direct potentials in the four differential equations are 
evaluated and some change is made in the limits of the exchange potential integrals, these equa- 
tions can be written in the following matrix form: 


_d^ 

dr 2 


+ k 2 - 


On + 1) 


“( k n 


0 


2V«(k n Z n , r ) . 


(23) 


The four components of u are the four radial functions of the free electron, v is a 4x4 symmetric 
matrix that is the sum of three matrices, 


•Reference 8, Equation 3.14. 




>- 




(24) 
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where is the direct potential, E.. is the exchange potential, and where both are functions of r . 
The matrix E. ^ contains in addition integrals with respect to r , and for the purpose of numerical 
integration it can be written as the sum of two matrices. The explicit forms of D. . , F ij? g* , 
and h.K are given in Appendix B. The value of a is 2 for i = j = 3 and i = j = 4, and is 1 for all 
other values of i and j . It is understood that for the exchange terms the components of « on the 
right-hand side of Equation 23 are inside the integrals of the exchange terms. 

Derivation of the Transmission Matrix from Solutions 
of the Differential Equations 

The method is similar to that used by Bransden and McKee (Reference 10), and by Marriot 
(Reference 1). Equation 23 constitutes a set of four coupled second-order differential equations. 
Three components of u can be eliminated from these equations, resulting in an 8 th order differ- 
ential equation for the remaining component. Therefore there are eight sets of solutions to Equa- 
tion 23. However, only half of these solutions are regular at the origin. Each of the four regular 
solutions corresponds to a definite vector u . The four vectors can properly be represented by a 
4x4 matrix u nj , n, j = 1,2,3, 4, where n corresponds to the particular component and j corre- 
ponds to the particular solution of u. The four solutions are carried out numerically in the next 
section. 

From the explicit form of V it can be seen that v vanishes at infinity. The asymptotic solu- 
tion of u as given by Equation 23 is therefore 

U nj(0 - a nj Sin ( k n r " ~Y~ + S nj ) <n, j = 1 , 2, 3, 4) , (25) 

where a nj is the amplitude and S nj is the phase shift of the j th solution of the n th component of u. 

Corresponding to the four components of u, there are 4 channels open to the reaction. If the 
incident wave is in the m th channel (m * 1, 2, 3, 4), the traveling wave in the n th channel will be 
given by 



/ 1 \" 


/ 1 , \' 


U n (r) % exp 

/ i ( k „ r ' T l n n ) 

S(m, n) - S mn exp 

H k n r_ 7 l n n ) 

(n = 1, 2, 3, 4) 


The constants S mn are the amplitudes of the scattered waves. Since Equation 26 is also the 
asymptotic solution of Equation 23, they must be equal to linear combinations of Equation 25. If 
we call the coefficients of the linear combinations P j , we must have 


4 



j=i 


r 

. /, 1 \1 

„ „ r. /. i 

1 

(exp 

' l ( k « r_ 7 l „ m )J 

5( m . n ) - ex P[ 1 l k " r ' T n 71 j 

j 


(n, m = 1, 2, 3, 4) . (27) 
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On the right-hand side, we have used the normalization of Blatt and Weisskopf (Reference 11). If 
we equate the coefficients of exp|^- i (k n r - 1/2 Z n ^)]and exp Ji (k n r - 1/2 l n in Equation 27, we 
obtain 


L 

j=i 


P a 

J nj 


exp - [i S nj ] - 


- 2i 


S(m. n) 


4 



j=l 


Separation of Equations 28 into real and imaginary parts gives 




> 


J 


(28) 


4 


i=i 




S(m, n) , 




4 

J~|[(RP i )co s 8 n ^(AP.) sinS nj ]a n . = 0, 

j=l 


4 

^j( RP ,) sillS nj + ( aP i) COsS nj] a nj 

j=l 



(29) 


4 


3=1 


cosS nj - (A p j) 


Sin8 nj 





J 


In the above R or represent the real or the imaginary part of the quantity that follows them. 
Equations 29 are a set of 16 linear equations for 16 unknowns RP jf AP jf Rs^ and £S mn . Once these 
unknowns are found,* the magnitude of S mn will be given by 


S I 2 

mn 


= (BsJ a +( As -J a 


( 30 ) 


* Equations 29 with their present form and without further simplifications are solvable by the computer. 
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The cross section is obtained by asymptotic expansion in spherical harmonics of the incident 
plane wave* 



The magnitude of the ingoing wave on the right-hand side of Equation 27 for n = m is 
[ k mA { 2l m + *)] 1/2 t* mes tJie magnitude of the partial wave of the expansion of r exp [ik m z ] . The 
plane wave has a flux of v which, in atomic units, is equal to k . The ingoing flux of the right-hand 
side of Equation 27 is therefore k m 2 /\jr[2l m + l)]. The outgoing flux in the channel n / m is |s mn | 2 . 
The cross section is obtained when we average the ratio of the outgoing flux to the ingoing flux 
over the initial states and sum over the final states. For a particular spin state of the two elec- 
trons, unpolarized electron beam, and unoriented atoms, the multiplicity of the initial states is 
(2 l t + l) (2 1 2 + l) , where l x and l 2 are the angular momentum of the bound and free electrons. 

For a polarized beam, m 2 = 0, where m 2 is the magnetic quantum number of the free electron. 
Then m 1 = M, where m x and M are the bound electron and the total magnetic quantum numbers. 
Since M is constant, there is only one initial state for a polarized beam. The multiplicity of the 
final states is 2 L + i, where L is the total orbital angular momentum. Since = l 2 , the cross 
section for m t n is 


Q 

'ran 


77 (fl, + l) 

k * 2 ( 2 Z l +1 ) 


IM 2 <« * n) 


(32) 


The outgoing partial wave in the incident channel m consists of the scattered wave plus the out- 
going wave given in the expansion of the plane wave. Then, according to Equation 27 for n = m , 
the magnitude of the amplitude of the scattered wave is 1 1 - S mn j . The elastic scattering cross 
section is therefore given by 


Q 

'mm 


77 ( 2 L + l) 

K 2 Wi 7 i) 


1 -s 2 

mm 


If we define a matrix T by the relation 


(33) 


T = 1 - S , 


(34) 


•Reference 1 1 , Ch. VIII, Equation 2.7. 
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Equations 32 and 33 can then be combined into a single equation, 


Q 

n 


77 ( 2L + 1 ) 

+ 1) 



(35) 


T mn is the transmitted amplitude in the n th channel due to an incident wave in the m th channel. 
The elements of T mn constitute the transmission matrix. 

The matrix S has two properties that are useful as tests on the accuracy of numerical inte- 
gration. From Equation 26 it can be seen that S transforms the ingoing wave into the outgoing 
waves. The continuity of the electronic current requires that S be a unitary matrix 


4 

22 K"i 2 = 1 (m = 2 ’ 3 ’ 4) ■ w 

n=l 


Furthermore, since the Hamiltonian is Hermetian, s must be symmetric (Reference 11): 


S 


mn 


s 


nia 


(37) 


Equations 36 and 37 are used as tests on the accuracy of numerical integration. 


A Useful Relation 


A relation based on the symmetry of the interaction potentials, which serves as another test 
on the accuracy of the solutions, can be derived. The l** 1 and the k th solutions of the i th com- 
ponent of u by Equation 23 are given by 


d 2 

Ldr 2 


+ k 2 - l. 





l. 

1 


(*i +1 )' 


Z V ‘i U i* ’ 

j 

j 




> 


j 


(38) 


Multiplying the first by u ik and the second by u 4l , subtracting the two expressions, and summing 
over i gives 



22 v ii [ u ik u iz- u ii u ik] • 


(39) 
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Since W i} = . , the interchange of the summation indices changes the sign on the right-hand side 

of the equation; the right-hand side must therefore be zero. Integrating the left-hand side from 
zero to infinity, we obtain 



(40) 


Integrating the above equation by parts, and applying Equation 25, we obtain 

4 

J - " kj a. k a ; j sin (S ik - 8^) = 0 (k, l = 1, 2, 3, 4, k t l) . (41) 

i=l 

Although the terms containing the exchange potentials do not cancel out on the right-hand side of 
Equation 39, the cancellation does take place after the integration is carried out in Equation 40. 


Transmission Matrix According to Born Approximation 

The Born approximation consists of neglecting the exchange potential terms appearing in the 
V matrix of Equation 23, and also of neglecting all the direct potential terms in this matrix except 
those terms that connect the incident channel to all other channels (Reference 12). Equation 23, 
when the incident wave is in the m th channel, reduces to 


dr 2 


+ k 




= 2D u 

n m m 


(n 


1, 2, 3, 4) 


(42) 


u m and u n are given asymptotically by 


U m - K W2 Sin ( k m r " ! .l) ’ (43) 

U n % k n ' 1/2 B nm cos (k n r - l„ 2 ) • (44) 

We have chosen the constants of proportionality of u m and u n such that B nm is the Born approxima- 
tion of the reactance matrix R (Reference 13).* Equation 43 shows that u m must have the follow- 
ing form (Reference 14): 


= k 


1/2 


( k * r ) 


(45) 


*Also Reference 11, Ch. X, Sec. 4. 
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where j x (k ffl r) are spherical Bessel functions. Furthermore, if y n represents the homogenous 
solution of Equation 42, it must have the following forms: 


Vn - a n k n r j !n (k„f) 


with a n some unknown constant. Multiplying Equation 42 on the left by y n and integrating the re- 
sult from zero to infinity, we obtain by partial integration 


2 y D nm u dr 

1 J n n m m 


V + k 2 

y ° dr 2 n 


K ( Z n + l) 


d d "I ' 

Vn d? U n “ U n dF y n 
U -J , 


The last equality has been obtained by noticing that y n and u n vanish at the origin, and by using 
their asymptotic forms as given by Equations 44 and 47. We therefore have 


2 ( k n k m ) 1/2 j Hn ( k n r ) D n m k* ( k m r ) r2 dr 


This is identical to the expression given for B by Seaton.* 

The transmission and the reactance matrices are related by T - - 2iR/( 1 - iR). Since in the 
Born approximation R = B « 1, the transmission matrix according to the Born approximation is 
given by 


T® = 4 

n m 


1 ( k „ k m) 1/2 j U» ( k n 0 D nm Urn (K r ) ^ dr 


Substitution of Equation 49 into 35 would give the cross section according to the Born approxima- 
tion. It should be noted that the symmetry of T insures Equation 37 to be satisfied while Equation 
36 is no longer satisfied. 


♦Reference 13, Equation 3-10. 
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Elastic Scattering of Electrons by the 2p States of the Hydrogen Atom 


The angular momentum of the free electron l 2 in the 2p channel has the values L - l , L , L + 1 , 
where L is the total angular momentum of the system. The first and the last values were consid- 
ered in previous sections. The case l 2 - L corresponds to a wave function in the 2p channel with 
a parity different from all channel wave functions considered previously. It therefore corresponds 
to elastic scattering. The wave function in this case is given by 


( F l* r 2) 



c 1LL 


^2p mi 


( r 0 


u(k 2 L, r 2 ) 


»i (° 2 ) 


(50) 


When Equation 9 is formed with this wave function, and minimized with respect to u(k 2 L , r 2 ) , 
treatment which led to the derivation of the four differential equations will give the following dif- 
ferential equation: 


+k 2_ k l + i> + 

.dr 2 2 r 2 


2 " 

Tj u l 


(r) 


+ k 2 2 ) S(L. 1 ) r Ejj (r) (2p|k 2 l) 


+ 2 


y 0 (2p2p, r) - jy 2 (2p2pi r) ] UL ( r) 


+ 2/3rR 21 (r) 


3 y L -i (2p k 2 L, r) 3y L+ , (2pk 2 L, r) "I 
( 2L + 1 ) ( 2L - 1 ) + (2L+ 1) (2L+ 3) J " 


(51) 


The asymptotic solution of this equation is given by 


u L % a L sin (k 2 r -L tt/2+ S L ) . ( 52 ) 

If the scattering amplitude is designated by T ss , it can be shown from the section on page 8 that 
for a particular L 


T 55 ~ 1 - exp 2i S - - 2i exp ( iS ) sin S , (53) 

where for simplicity we have suppressed the subscript L. The cross section, according to Equa- 
tion 35, is given by 


4? t ( 2L +1 ) 


3k 2 2 


(54) 
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The total elastic scattering cross section by the 2p states is the sum of this cross section and 
the cross section corresponding to i 2 ~ l- 1 and l 2 = L + l (which were considered previously). 

The Born amplitude, Equation 48, in this case is given by 


B 


55 



where, by Equation 51, 


+ y 0 (2p2p, r) (2p2p, r) 


(55) 


(56) 


NUMERICAL INTERGRATION 

Decomposition of the Differential Equations 

If it were not for the definite integrals appearing in the potential matrix v , the set of the four 
coupled differential equations (23) could be integrated by any standard technique. The presence 
of these unknown constants whose integrand involves the unknown functions makes it necessary 
to solve these equations by iteration or by transformation of M into other vectors, whose differ- 
ential equations do not contain definite integrals. Since the terms containing definite integrals 
are small as compared with the direct potentials, the iteration method can be used by assuming 
that the values of these integrals are zero. The differential equations are then integrated, the 
values of the definite integrals that are subsequently obtained are substituted in the differential 
equations, and the integration is repeated. The process is repeated until sufficiently consistent 
values of these integrals are obtained. This method is useful if the convergences of the constants 
are fast enough, and the cross section is not very sensitive to the values of these constants. 

In the second method, the transformation of u fixes the values of the constants and thus 
avoids iteration, whereby the computation is reduced considerably. The description of the method 
will be given here (Reference 15).* 

By making use of Equations 24, Equation 23 can be written 

4 


j = l 




L 

M=1 


g/t 


C *1 


(57) 


+ k > 
dr 2 ‘ 


l i i l i +1 ) 


•Also, Reference 1. This description differs from the description of Reference 15 and the present paper. In Reference 1, Vj in Equa- 
tion 61 is set to 0; this makes - 0. Equations 62 then reduce to a set of homogeneous equations whose determinant must be 0. 
Since the amplitude of any of the four components ofu can be left arbitrary, one of the is set equal to 1 and the rest of the con- 
stants are found subsequently. 
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where 


C/f 

i j 


J o 


(r)dr . 


(58) 


We introduce the functions v. and u> 1 that are solutions of the following differential equations: 

4 

= 2 2Z [ D d +F ij] v i • (59) 


d 2 , „ 
^ +k i 2 - 


dr 2 ‘ r 2 


, kl 


2 +Fi 'i u j k! + 28 < i - k > 

i=i 


*kl 


(60) 


Then u i is given by the following expression: 


4 4 

u * = v i + XI H Z ] C * 1 u * k! • (61) 

k=l 1 = 1 V=l 

Equation 61 can be verified by multiplying Equation 60 by C^; summing over k, i, and v; and add- 
ing to Equation 59— whereupon Equation 57 results. Substitution of Equation 61 in Equation 58 
gives 

4 4 cr 

L H H [s ( i)M,kZv)-A/f]c^ = BM 

k=l 1=1 V=1 

(i,j - 1,2,3, 4; \± - l, 2 for i = j - 3 and i = j ~ 4; - l otherwise) > (62) 

where k^ 1 andB.^. are defined by 

ij i j j 

A /* 1 = h/* u kJ dr , 

‘ J Jo “ ’ 

B i1 = f h t s ^ 

Jo 

The numerical integration is carried out by integrating Equations 59 and 60 by any standard 
method, calculating k{* 1 and B.^ by Equations 63 and, finally, solving the system of 18 algebraic 
equations given by Equation 62 to find C * t . With the known values of these constants the integra- 
tion of Equation 57 is straightforward. 


(63) 
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The determinant of Equation 62 becomes singular for L - 0 and 1. This is shown in Appendix 
C. To remove the singularity, some of the C k ^ are chosen arbitrarily, and the rest of the are 
found in terms of the chosen ones. 

Solution at the Origin 

In order that the four solutions of u be independent of each other, we must have 

4 

2] C i % * 0 (^1-2, 3, 4) , (64) 


where c. are some constants. A necessary condition for this to be satisfied is that the deter- 
minant of Equation 64 be nonzero: 


I KII * 0 • (65) 

It is not difficult to see that this also is a sufficient condition. At the origin the solution u.. can 
be expressed as power series in r , 


CD 

■ L 


( 66 ) 


where a.*j are the coefficients of expansion, and s. are given integers for each component of u 
and are fixed by the behavior of Equation 23 at the origin. We can satisfy Equation 65 at the origin 
by having 


IKII * 0 . (67) 

By choosing suitable values of a. 0 . , subject to the restriction of Equation 67, four independent 
solutions are obtained. 

Solution at Large r 

With given initial values, the solution of Equation 23 can be extended from origin to any de- 
sired value of r . To obtain the asymptotic amplitudes and the phase shifts, the presence of the 
centrifugal and long-range potentials makes it necessary to extend the solutions to infinity. This 
is undesirable because of the time consumption on the computer, and the accumulated error due 
to the long-range integration. Seaton (Reference 13) has solved the problem of r" 2 long-range 
potentials occurring in the off-diagonal terms of the potential matrix v by diagonalizing the 
asymptotic form of the differential Equation 23 and the corresponding s- matrix. By an inverse 
transformation the elements of the original S -matrix are found. 
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Instead, we develop here a perturbation theory which is based on the method described by 
Mott and Massey (Reference 16). The error in the resulting solution is inversely proportional to 
the square of the distance from the origin. 

Equation 23 for large distances of r can be written 


_di 

dr 2 



0 


2 U„(k n i n ,r) . 


( 68 ) 


where u is the sum of the centrifugal potential matrix and the asymptotic form of the v matrix. 

The elements of u are given in Appendix D. A component of Equation 68 is of the following form: 

^2 +k2 ] u ( r ) ^ e( r ) ’ 

g(r) « k 2 u(r) (g(r) - 0 as r -» oo) 

The perturbation theory is applied between some large distance R and infinity. Suppose u vanishes 
atR; then we have the following boundary condition: 


(69) 


u(R) = 0 . 


(70) 


If we represent the solution of the homogenous equation by y(r), at infinity, we must have 


y(r) - asin(kr-kR) , 

u(r) = (a + Aa) sin (kr - kR + 77 ) 


(71) 


where a is the amplitude of u(r) if g(r) were identically zero and AA and 17 are generated by 
g(r). Since g(r) is small, we can write 

u = y(l + O • (72) 


where ( is a small function. Substitution of Equation 72 in Equation 69 gives 


_d_ 

dr 



g(r)y . 


(73) 


where, upon double integration, we obtain 


C 



r 


g(r')ydr' 


( 74 ) 


The constants of integrations are fixed by the condition (70) and the fact that u' (R) = y' (R) . 
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We now integrate Equation 74 by parts, 


f g( r ) y dr f ^ 
J R J R y 


f g(r)ydr f * 

j r j r y 


(75) 


When the integration with respect to y is carried out, and the result is substituted in Equation 73, 
we obtain 


u(r ) - sin (kr - kR) 


a + 


1 f r 

IT J g(r) cos (kr - kR) dr 

^ D 


+ cos (kr - kR) 




g( r ) sin (kr - kR) dr 


Comparison of the second of Equations 71 and 76 shows that 


Aa 


_ac 

= k" J g(r) cos (kr - kR) dr , 


1 f °° 

V ~ “ afc J g(r) sin (kr - kR) dr 

R 

to first order. The functions g(r) in the four differential Equations 68 are given by 

MO = 2 22 b ijUj . 

j 

To first order this can be written by 

MO = 2^]a j U ij sin(k j r-k jRj ) , 


(76) 


(77) 


(78) 


(79) 


where R. is the last zero of vr with positive slope. Substitution of this equation in Equation 77 
gives 


Aa . 


- T.T.j cos ( k i r ~ u i R i) U ij sin ( k i r-k. Rj ) dr , 
) 

J2 | ‘ sin (k ; r - k; Rj) U ij sin (k ; r - Rj Rj ) dr 


( 80 ) 
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A a . and rj. can easily be calculated by substituting the values of u.. from Appendix D, integrating 
the resulting integrals by parts, and retaining the leading terms. 

The asymptotic amplitudes and phase shifts are given by 


a. (oo) - a . (R.) + A 9i , 

S i (°°) r s i ( R i) + t?. + [L- S(i, 3) + S(i, 4)] 3 * 


(81) 


where a. (r.) and 8 . (R i j are the amplitudes and total phase shifts calculated at r. by the ma- 
chine, and where S(i , 3) and 8 (i , 4) are the 8 functions. 


Details of the Numerical Integration 

Milne's (Reference 17) method with variable mesh size and Simpson's* rule were used for 
the integration of the differential equations and evaluations of the integrals respectively. As the 
solution advances from the origin, the differential equations become less sensitive to the size of 
the increment, and the error of integration falls below certain small number e . At each value of 
r the value of the function is found, first with the given value of the increment, and second with 
the value of increment divided in half. The error of integration is defined as the difference be- 
tween these two solutions. When the error becomes small, the increment is doubled until a 
maximum value is reached. At some distance R t all the exchange potentials and, similarly, all 
the direct potentials except those representing optically allowed transitions and the 2p - 2p elastic 
scattering potential become vanishingly small (see Appendix D). At this distance the set of dif- 
ferential equations is replaced by the simpler set containing only these potentials. The integra- 
tion is continued until some distance R 2 , where the first-order solution of the rest of the range of 
integration is obtained by the method developed in the previous section. No attempt was made to 
solve any set of linear equations or any matrix equations, as these equations are solvable by the 
computer in their original form. 

The values of the constants of the numerical integration are given below; h. and h f are the 
initial and the final increments of integration. In some exceptional cases, different values were 
used. 

h A = l x 10~ 5 

h f = 0.05 

e = 1 x 10" 4 
R x = 30 

R 2 = 200 

(All quantities are in units of Bohr radius except e, which is dimensionless.) 


*Reference 17, Sec. 33. 
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RESULTS AND DISCUSSION 


The differential equations (23) with the known values of the elements of the potential matrix 
V as given in Appendix B were integrated numerically by the methods described in the last sec- 
tion. By choosing different values for the determinant (Equation 67) different sets of independent 
solutions can be generated. The cross sections reported in this paper have been obtained by 
averaging the cross sections obtained from two independent sets of solutions. To test the ac- 
curacy of the numerical integration, we define the three quantities D mn , D m ' n , and D m given by 




( 8 . - B. ) 

\ im in / 


* 

/ k a . a sin (8 . “8.1 

/ i im m \ un in/ 


(m, n - 1, 2, 3, 4, m / n) , 


(82) 


D ' 


m n 


s - s 

i mn nm | 


Is i 

+ 

Is i 

1 mn | 


1 nm 


(m, n - 1, 2, 3, 4, m / n) 


D 

m 



(m = 1, 2, 3, 4) . 


(83) 


(84) 


Based on Equations 36, 37, and 41 in an exact solution of the four differential equations, the right- 
hand side of these equations would vanish; they can therefore be used to test the accuracy of the 
numerical integration. As an illustration the numerical values of D , D ' , and P for the case 

mn 7 mn 7 m 


2s - 2p coupling, /S = + 1 , k 1 

= 2.0, and L = 

3, are given below: 



D . 2 

= i.4 x i<r 3 

D !3 = 

2.6 x 10' 4 , D 14 = 

1 . 3 x 

10" 3 , 

D 23 

= 5.1 X 1(T 4 

d 24 = 

2.2 x 1CT 3 , d 34 = 

1.8 x 

10~ 3 , 

°12 

= 7.6 x 10’ 4 

D ,3 = 

5. 1 x 1(T 3 , D 14 = 

5.6 x 

10‘ 3 . 

D 23 

a 5.4 x 10~ 3 

D 2'4 = 

5.7 x 10' 3 , D 3 ' 4 = 

1.3 x 

10' 3 , 

D t " = 1.8 

x 1(T 4 , D 2 " 

= 7.6 x 10' 5 

, D/ = 2.5 x 10’ 4 , 

d; ^ 

= 4.4 X 


To compare the results of the numerical integration by the noniterative method we have 
carried out here with those of the iterative method of References 2, 4, and 5, we have provided 
Table 1.* The is - 2s excitation cross section is given by the two methods; I and II refer to the 


•Author is indebted to Dr. K. Smith for sending some of the data used in Table 1. 
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Table 1 


iterative and noniterative methods, 


Comparison of the Iterative and the Noniterative Results 
for the Singlet (L = 0, 1; k 4 = 0.9, 1.0), Is- 2s Excitation 
Cross Section.* 


(a) Is - 2s coupling 


k ! 

L 

- 2 , 

E „ a „ (percent) 

n 

II 

I 

II 

0.90 

0 

0.0384 

0.0375 

7.1 

0.72 

1.00 

0 

0.0714 

0.0725 

unknown 

0.53 

0.90 

1 

0.008 

0.0017 

386 

0.91 

1.00 

1 

0.051 

0.0583 

55 

0.75 


(b) Is - 2s - 2p coupling 


k i 

L 

Qh- 

- 2s 

E m,x (percent) 

*1 

II 

I 

II 



0.90 

0 

0.0529 

0.523 

0.40 

0.40 

1.00 

o 

0.0766 

0.0768 

0.12 

0.60 

0.90 

1 

0.0045 

0.0048 

2.3 

10 

1.00 

1 

0.0145 

0.0147 

0.33 : 

1.3 


*1 and II refer to iterative and noniterative methods, respectively; E max is the 
maximum of the error to value ratios in the reciprocity relation. 
tSee Reference 2. 
tSee References 4 and 5. 


respectively; and E max is the maxi- 
mum of the error to value ratios in 
the reciprocity relations (Equation 
83). In the is - 2s eigenstates 
coupling approximation the non- 
iterative method is far more ac- 
curate than the iterative method 
and, as is seen, the cross sections 
by the two methods differ from each 
other sometimes in their first 
significant figure. In the is - 2s - 2p 
eigenstates coupling approximation, 
on the other hand, the results by 
the iterative method seems to be 
somewhat more accurate. The 
reason is contributed to the effect 
of the r" 2 long-range potential, 
which appears in the differential 
equations when the 2p state is in- 
cluded in the eigenstates coupling 


approximation. Two different methods are used in References 4 and 5 and in the present paper 


to estimate the effect of this potential for large distances; and it may be that in References 4 and 


5 this effect is better taken into account. Nevertheless the cross sections are the same in their 


first three decimal places. 


In Figure 1 we present the theoretical and the experimental estimate of the is - 2s excitation 
cross section. The calculated curves are Born, is - 2s coupling, is - 2s - 2p coupling exchange 
neglected, and Is - 2s - 2p coupling exchange included approximation. The first three of these 
curves are the same as References 4 and 5. The experimental curves are those of Lichten and 
Schultz (Reference 18) and Stebbings, Fite, and Hummer (Reference 19). The various calculated 
results agree better with the results of Lichten and Schultz. However, recent calculations of 
Taylor and Burke (Reference 20) have shown that, in an eigenstates expansion calculation where 
Is, 2s , 2p, 3s , and 3p are included, the cross section at the peak of the is - 2s - 2p curve is re- 
duced by 30 percent. This suggests that, within the eigenstates expansion approximation, more 
states should be included to insure that the convergence has been achieved; and the discrepancy 
between the two experimental results is still an unresolved problem. As another theoretical ap- 
proach to the problem, H. L. Kyle and A. Temkin (Reference 21) have extended the nonadiabatic 
theory of scattering developed by A. Temkin (Reference 22) to the L = 0 , Is - 2s inelastic scatter- 
ing of electrons by the hydrogen atom. They find a 30 percent decrease in the Is - 2s cross section 
as calculated by the is - 2s close coupling approximation. 

Comparison of the exchange neglected and exchange included is - 2s - 2p coupling shows that 
exchange is mostly important at threshold, and its effect does not extend beyond 20 electron volts. 
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Table 2 gives the numerical values of 
the is - 2s cross section in different approxi- 
mations. The Is - * 2s excitation cross section 
in the singlet state has an interesting be- 
havior immediately above threshold. In 
Figure 2 this cross section for a range of 
600 milli electron volts (mev) above 
threshold is plotted. In the is - 2s coupling 
approximation a maximum appears at 34 
mev, while in the Is - 2s - 2p coupling ap- 
proximation there are three maxima of 
approximately the same magnitudes at 17, 

34, and 87 mev, respectively. In the singlet 
case the cross section rises sharply within 
a range of 17 mev above threshold to a value 
of about 0.0477 a 0 2 ; it then rises with an ap- 
proximately constant and small slope. The 
contribution of the triplet case is seen to be 
almost negligible at the threshold, and it has 
no maximum in this region. It should be 
noted that the principal maximum in the 
is - 2s excitation cross section appears at 
about 3 ev with a value of about 0.35, and has 
a contribution from a higher angular mo- 
mentum than L = 0 . Although no study has 
been made to relate the existence of the maxima above threshold to any physical phenomena, it 
may be said that, similar to resonances below threshold in the elastic scattering of electrons by 
the hydrogen atom, these maxima are due to formation of some unstable states of the negative 
hydrogen ion. The numerical values of the is - 2s cross section at threshold are given in Table 3. 
Damburg and Peterkop (Reference 6), and Gailitis and Damburg (Reference 23) have made an 
extensive study of the behavior of different cross sections near threshold in the is - 2s and the 
Is - 2s - 2p eigenstates coupling approximations. 

In Figure 3 we have shown the 2s-* 2s elastic cross section. The is - 2s coupling approxima- 
tion gives a value of 944 77 a 0 2 at zero incident energy, while the corresponding value in the Born 
approximation is 768 77 a 2 . The high value of this cross section at zero energy is in sharp con- 
trast with its geometrical cross section. The zero energy 2s - 2s cross section in the is - 2s - 2p 
coupling approximation, because of the r 2 potential, is difficult to find. The 2s - 2s cross section 
has certain maxima and minima at low energy which are not found in the is - is cross section. 
Figure 4 shows the L = 0 singlet and triplet 2s - 2s cross section in the two approximations. While 
there is one minimum in the is - 2s coupling approximation, there are three minima in the 
is - 2s - 2p coupling approximation. It is thought that the existence of these minima is due to a 
wider potential range in the 2s -2s scattering, a case which does not exist in the is- is 



Figure 1 — Is - 2s total excitation cross section. (Is - 2s 
refers to Is - 2s eigenstates coupling approximation; 
Is - 2s - 2p has similar meaning. EX. NEGL. refers to 
exchange neglected case, BORN is the Born approxima- 
tion, and EXP. refers to experiment.) 
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Table 2 

Is - 2s Excitation Cross Section, 


(a) Born approximation 


ki 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L - 7 

2 

Qt 

0.9 

0.16376 

0.00981 

0.00021 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.17378 

0.17379 

1.0 

0.19578 

0.04795 

0.00428 

0.00026 

0.00000 

0.00000 

0.00000 

0.00000 

0.24827 

0.24827 

1.1 

0.16272 

0.07073 

0.01141 

0.00125 

0.00011 

0.00001 

0.00000 

0.00000 

0.24622 

0.24623 

1.2 

0.12704 

0.07896 

0.01858 

0.00299 

0.00039 

0.00004 

0.00000 

0.00000 

0.22800 

0.22800 

1.5 

0.05872 

0.06606 

0.02979 

0.00939 

0.00242 

0.00053 

0.00010 

0.00002 

0.16703 

0.16706 

2.0 

0.01946 

0.03363 

0.02521 

0.01365 

0.00614 

0.00236 

0.00081 

0.00025 

0.10151 

0.10187 

3.0 

0.00388 

0.00909 

0.01019 

0.00866 

0.00628 

0.00394 

0.00226 

0.00116 

0.04546 

0.04758 

4.0 

0.00123 

0.00320 

0.00420 

0.00431 

0.00385 

0.00301 

0.00216 

0.00141 

0.02337 

0.02720 


(b) Exchange neglected Is - 2s - 2p eigenstates coupling approximation 


ki 

o 

II 

L = 1 

L = 2 

L = 3 

II 

-1 

L = 5 

L = 6 

n 

2 

Qt 

0.9 

0.2202 

0.0749 

0.3535 

— 

-- 

— 

— 

— 

0.6486 

0.6486 

1.0 

0.1685 

0.1427 

0.1598 

0.0517 

— 

— 

— 

— 

0.5227 

0.5227 

1.1 

0.0951 

0.1142 

0.0298 

0.0616 

0.0231 

— 

— 

— 

0.3238 

0.3238 

1.2 

0.0594 

0.1137 

0.0032 

0.0360 

0.0244 

0.0135 

— 

— 

0.2502 

0.2502 

1.5 

0.0249 

0.0861 

0.0201 

0.0068 

0.0112 

0.0118 

0.0074 

— 

0.1683 

0.1683 

2.0 

0.0101 

0.0373 

0.0255 

0.0107 

0.0046 

0.0034 

0.0033 

— 

0.0949 

0.0953 


(c) Is - 2s eigenstates coupling approximation 


Singlet 

ki 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2 S 

— 

0.9 

1.0 

1.1 

1.2 

1.5 

2.0 

0.0375 

0.0725 

0.0701 

0.0547 

0.0241 

0.0072 

0.0017 

0.0583 

0.0525 

0.0534 

0.0384 

0.0157 

0.0000 

0.0002 

0.0023 

0.0054 

0.0110 

0.0093 

0.0000 

0.0000 

0.0000 

0.0002 

0.0022 

0.0041 

0.0000 

0.0000 

0.0000 

0.0000 

0.0004 

0.0015 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

0.0005 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0002 

— 

0.0392 

0.1310 

0.1249 

0.1137 

0.0762 

0.0385 


Triplet 

2 s + 

Qt 

ki 

L = 0 

L = 1 

L - 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 


0.9 

1.0 

1.1 

1.2 

1.5 

2.0 

0.0004 

0.0021 

0.0044 

0.0061 

0.0073 

0.0049 

0.1686 

0.1528 

0.1052 

0.0737 

0.0355 

0.0162 

0.0060 

0.0446 

0.0568 

0.0576 

0.0406 

0.0205 

0.0000 

0.0021 

0.0068 

0.0114 

0.0174 

0.0143 

0.0000 

0.0001 

0.0005 

0.0015 

0.0050 

0.0074 

0.0000 

0.0000 

0.0000 

0.0002 

0.0012 

0.0032 

0.0000 

0.0000 

0.0000 

0.0000 

0.0002 

0.0012 

— 

0.1750 

0.2017 

0.1737 

0.1505 

0.1072 

0.0677 

0.2142 

0.3327 

0.2987 

0.2642 

0.1833 

0.1062 

0.2142 

0.3327 

0.2987 

0.2642 

0.1833 

0.1068 


(d) Is - 2s - 2p eigenstates coupling approximation 
Singlet 


k, 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

r* 

ii 

CD 

L - 6 

L = 7 


— 

0.9 

0.0523 

0.0048 

0.0620 

— 

— 

— 

— 

— 

0.1191 

— 

1.0 

0.0768 

0.0147 

0.0833 

0.0092 

— 

— 

— 

— 

0.1840 


1.1 

0.0585 

0.0245 

0.0647 

0.0236 

0.0055 

— 

— 

— 

0.1768 


1.2 

0.0382 

0.0251 

0.0246 

0.0252 

0.0081 

0.0028 

— 

— 

0.1240 


1.5 

0.0123 

0.0308 

0.0015 

0.0041 

0.0051 

0.0034 

0.0023 

0.0026 

0.0621 


2.0 

0.0049 

0.0152 

0.0068 

0.0021 

0.0010 

0.0010 

0.0011 

0.0008 

0.0329 


3.0 

0.0010 

0.0031 

0.0031 

0.0023 

0.0015 

0.0009 

0.0005 

0.0006 

0.0130 


4.0 

0.0003 

0.0010 

0.0012 

0.0012 

0.0010 

0.0008 

— 

— 

0.0055 



Triplet | 

2 S + 2t 

Q 

^ r 

ki 

L = 0 

L = 1 

L = 2 

CO 

II 

ii 

L = 5 

L = 6 

L = 7 

2 t 

0.9 

0.0013 

0.0748 

0.0019 

— 

— 

— 

— 

— 

0.0780 

0.1971 

0.1971 

1.0 

0.0040 

0.1224 

0.0195 

0.0214 

— 

— 

— 

— 

0.1673 

0.3513 

0.3513 

1.1 

0.0050 

0.1013 

0.0326 

0.0077 

0.0131 

— 

— 

— 

0.1597 

0.3366 

0.3366 

1.2 

0.0055 

0.0724 

0.0359 

0.0036 

0.0105 

0.0076 

— 

— 

0.1355 

0.2596 

0.2596 

1.5 

0.0045 

0.0333 

0.0309 

0.0072 

0.0046 

0.0054 

0.0049 

0.0043 

0.0951 

0.1573 

0.1573 

2.0 

0.0031 

0.0155 

0.0176 

0.0101 

0.0044 

0.0025 

0.0023 

0.0019 

0.0574 

0.0903 

0.0907 

3.0 

0.0013 

0.0048 

0.0065 

0.0059 

0.0044 

0.0029 

0.0018 

0.0012 

0.0288 

0.0418 

0.0439 

4.0 

0.0006 

0.0019 

0.0028 

0.0030 

0.0027 

0.0022 

-- 

— 

0.0132 

0.0187 

0.0261 
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CROSS SECTION (ir ag) 



Figure 2 — L = 0, Is - 2s excitation cross section above threshold, (The cross sections are given for the two spin states 
singlet and triplet, and for the two approximations Is - 2s and Is - 2s - 2p. The total cross section is the sum of the 
singlet and the triplet cross sections.)* 


scattering. In Table 4 we have listed the nu- 
merical values of the 2s — 2s cross section in 
different approximations. 

In Figure 5 the four calculated curves for 
the is — 2p excitation cross section are com- 
pared with the measurement of Fite, Stebbings, ^ 

CMC 

and Brackmann (References 24 and 25). The £ 

is - 2s - 2p and the Born curves are the same as z 

o 

in References 4 and 5, but the ls-2s-2p ex- £ 

change neglected and the is - 2p curves are not ^ 

calculated in these references. As concluded <§ 
before, the calculated curves are higher than 
the experimental. Moreover, we notice that, 
similar to the is-* 2s excitation cross section, 
the inclusion of the exchange lowers the value 
of the cross section at threshold. Table 5 gives 
the numerical values of the is-2p cross sec- 
tion in different approximations. 

The calculation of the 2p - 2p elastic 
cross section is more complicated than the 



Figure 3 — 2s -2s total elastic cross section. (Curves 
are designated as in Figure 1.) 


According to Gailitis and Damburg, when the energy difference between the 2s and 2p states is neglected in the ls-2s-2p couplings, 
the Is — 2s excitation cross section does not go to zero at threshold (see Reference 23, and Figure 2). Figures 2 and 4 show that in 
the Is - 2s - 2p couplings, if E n _i and E n represent the energy with respect to the threshold of the two neighboring maxima or minima 
then E n /E n -j = constant. This maybe attributed to the r~ 2 potential, which is due to the coupling between the 2s and the 2p states. 
For further details, see Reference 23- 
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Table 3 


The Singlet L = 0, Is - 2s Excitation Cross Section Near Threshold.* 


K 

0 

0.01 

0.02 

0.025 

0.030 

0.035 

0.04 

E(mev) 

0 

1.36 

5.44 

8.50 

12.2 

16.7 

21.8 

Q, 

0 

0.0168 

0.0298 

— 

0.0377 

— 

0.0423 

Qa 

0 

— 

0.0149 

0.0259 

0.0349 

0.0405 

0.0353 

k 2 

0.045 

0.050 

1 

0.060 

0.070 

0.080 

0.090 

0.100 

E(mev) 

27.5 

34.0 

49.0 

66.6 

87.0 

110. 

136. 

Q, 

— 

0.0446 

0.0441 

0.0435 

0.0423 

0.0412 

0.0405 

Qa 

0.0346 

0.0405 

0.0391 

0.0361 

0.0395 

0.0392 

0.0385 


*k2 is the wave number of the inelastically scattered wave, and E is the corresponding energy in mev; Qj and Q 2 are the cross sec- 
tions according to the Is - 2s and the Is - 2s - 2p couplings, respectively. 
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Table 4 

2s - 2s Elastic Cross Section. 


(a) Born approximation 


k ? 

L = 0 

L = 1 

L - 2 

L = 3 

L = 4 

L= 5 

L = 6 

L = 7 

2 

Qi 

0.24 

389.97 

15.910 

0.27680 

0.00319 

0.00003 

0.00000 

0.00000 

0.00000 

406.16 

406.17 

0.50 

105.52 

32.129 

5.0152 

0.55306 

0.04846 

0.00347 

0.00022 

0.00001 

143.269 

143.276 

0.68 

47.381 

24.746 

7.4981 

1.7023 

0.31602 

0.04869 

0.00658 

0.00078 

81.700 

81.703 

0.83 

26.562 

17.895 

7.5601 

2.5080 

0.69820 

0.16365 

0.03393 

0.00621 

55.427 

55.440 

1.23 

8.2461 

7.6846 

4.9702 

2.7483 

1.3500 

0.57745 

0.22452 

0.07803 

25.879 

25.990 

1.80 

2.4331 

2.7773 

2.2763 

1.6805 

1.1541 

0.71994 

0.41776 

0.22120 

11.680 

12.105 

2.87 

0.5245 

0.71381 

0.68903 

0.61123 

0.51876 

0.41048 

0.31076 

0.21860 

4.0011 

4.8280 

3.91 

0.18794 

0.27611 

0.28658 

0.27339 

0.25098 

0.21580 

0.17946 

0.14029 

1.8106 

2.7417 



(b) Exchange neglected Is 

- 

2s - 2p eigenstates coupling approximation 


k 2 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2 

_Sx 

0.24 

32.02 

8.489 

42.26 

— 

— 

— 

— 

— 

82.77 

82.78 

0.50 

2.041 

7.710 

23.35 

9.004 

9.368 

7.702 

5.870 

— 

65.045 

65.052 

0.68 

2.255 

8.247 

15.16 

5.547 

4.373 

3.755 

2.975 

— 

42.312 

42.316 

0.83 

1.861 

8.206 

10.67 

4.637 

2.792 

2.223 

1.818 

— 

32.207 

32.226 

1.23 

1.716 

4.987 

5.030 

3.321 

1.797 

1.058 

0.7179 

— 

18.537 

18.726 

1.80 

1.020 

2.134 

2.150 

1.750 

1.275 

0.8320 

0.5457 

— 

9.707 

10.352 


(c) Is - 2s eigenstates coupling approximation 


Singlet 


K 

L = 0 

L = 1 

L = 2 

CO 

II 

r 

ii 

4* 

L = 5 

L = 6 

ii 

2 S 

— 

0.24 

0.3303 

8.196 

0.2628 

0.0028 

— 

— 

— 

— 

8.792 

— 

0.50 

1.532 

10.38 

0.0275 

0.0048 

0.0008 

0.0002 

0.0249 

— 

11.97 


0.68 

1.115 

5.536 

1.502 

0.1150 

0.0087 

0.0010 

0.0017 

— 

8.279 


0.83 

0.8980 

3.512 

1.997 

0.4303 

0.0747 

0.0129 

0.0032 

— 

6.928 


1.23 

0.5702 

1.413 

1.236 

0.7010 

0.3129 

0.1228 

0.0450 

— 

4.401 


1.80 

0.2825 

0.5370 

0.5285 

0.4193 

0.2931 

0.1863 

0.1110 

— 

2.358 



L= 0 


Triplet 


L - 1 


L = 2 


L = 3 


L = 4 


L = 5 


L = 6 


L= 7 


2 t 




Q-r 


0.24 

0.50 

0.68 

0.83 

1.23 

1.80 


45.94 

0.2102 

1.366 

2.112 

1.811 

0.8989 


118.8 

34.44 

18.13 

11.65 

4.691 

1.735 


7.713 

21.05 

12.74 

8.725 

4.008 

1.652 


0.0540 

2.776 

4.059 

3.887 

2.585 

1.316 


0.2521 

0.8282 

1.230 

1.399 

0.9510 


0.0316 

0.1463 

0.3225 

0.6637 

0.6359 


0.0994 

0.0442 

0.0850 

0.2862 

0.4008 


172.51 

58.86 

37.31 

28.01 

15.44 

7.590 


181.30 

70.83 
45.60 
34.94 

19.84 
9.947 


181.31 

70.84 

45.60 

34.96 

20.03 

10.592 


(d) Is - 2s - 2p eigenstates coupling approximation 


Singlet 

k 7 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 


— 

0.24 

0.50 

0.68 

0.83 

1.23 

1.80 

2.87 

3.91 

7.800 

0.2858 

0.0661 

0.1675 

0.3739 

0.2416 

0.0847 

0.0365 

14.79 

0.6960 

1.044 

1.088 

1.075 

0.4974 

0.1489 

0.0635 

22.42 

3.149 

2.455 

2.105 

1.196 

0.5133 

0.1574 

0.0702 

4.447 

1.884 

0.9905 

0.7196 

0.4133 

0.1466 

0.0698 

2.925 

1.480 

0.8282 

0.3940 

0.2923 

0.1280 

0.0660 

2.063 

1.098 

0.6606 

0.2496 

0.1928 

0.1057 

0.0599 

1.491 

0.7928 

0.5071 

0.1827 

0.1247 

0.0852 

0.1409 

0.0799 

0.0658 

45.01 

15.057 

8.820 

6.347 

4.332 

2.3553 

0.9223 

0.3659 


Triplet 

S s + s T 


k 2 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

St 

0.24 

0.50 

0.68 

0.83 

1.23 

1.80 

2.87 

3.91 

16.52 

6.172 

2.346 

1.709 

1.391 

0.7898 

0.2622 

0.1113 

1.236 

17.19 

12.88 

9.166 

4.199 

1.656 

0.4652 

0.1945 

89.65 

30.69 

14.52 

8.976 

3.880 

1.603 

0.4857 

0.2146 

12.90 

7.937 

5.367 

2.681 

1.291 

0.4479 

0.2122 

5.559 

3.046 

2.365 

1.585 

0.9466 

0.3901 

0.2000 

5.373 

2.459 

1.514 

0.9110 

0.6520 

0.3234 

0.1812 

4.257 

2.090 

1.246 

0.5658 

0.4305 

0.2623 

0.4094 

0.2786 

0.2045 

107.41 

82.14 

45.28 

30.343 

15.622 

7.648 

2.8413 

1.1138 

152.42 

97.206 

54.099 

36.686 

19.955 

10.002 

3.7636 

1.4797 

152.43 

97.213 

54.103 

36.705 

20.066 

10.427 

4.5905 

2.7306 
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Table 5 

Is -» 2p Excitation Cross Section. 


(a) Born approximation 


ki 

o 

II 

L - 1 

L = 2 

L = 3 

L = 4 

r 

ii 

cn 

L = 6 

ii 

2 

Qi 

0.9 

0.00107 

0.46700 

0.09607 

0.01025 

0.00087 

0.00007 

0.00001 

0.00000 

0.56534 

0.57535 

1.0 

0.00499 

0.48867 

0.35645 

0.13467 

0.03964 

0.01032 

0.00254 

0.00059 

1.03787 

1.03851 

1.1 

0.00702 

0.36207 

0.41698 

0.24992 

0.11533 

0.04649 

0.01753 

0.00624 

1.22158 

1.22859 

1.2 

0.00747 

0.25540 

0.38137 

0.29903 

0.17908 

0.09269 

0.04453 

0.02004 

1.27961 

1.30741 

1.5 

0.00550 

0.09184 

0.19986 

0.23775 

0.21386 

0.16292 

0.11291 

0.07232 

1.09696 

1.28101 

2.0 

0.00234 

0.02222 

0.05938 ' 

0.09386 

0.11287 

0.11376 

0.10248 

0.08374 

0.59065 

1.04055 

3.0 

0.00048 

0.00285 

0.00806 

0.01537 

0.02309 

0.02920 

0.03271 

0.03273 

0.14449 

0.66256 

4.0 

0.00013 

0.00066 

0.00179 | 

0.00359 

0.00586 

0.00816 

0.01009 

0.01114 

0.04142 

0.45252 


(b) Exchange neglected Is - 2s - 2p eigenstates coupling approximation 


k, 

o 

ii 


L = 2 

L = 3 

L = 4 

L = 5 

CD 

II 

i'- 

ll 

*4 

2 

Qt 

0.9 

0.1600 

0.3985 

0.6497 

— 

— 

— 

— 

— 

1.2082 

1.2194 

1.0 

0.1007 

0.2917 

0.8190 

0.2190 

0.0476 

0.0224 

— 

— 

1.5004 

1.5041 

1.1 

0.0980 

0.2008 

0.6201 

0.3696 

0.1380 

0.0586 

0.0166 

— 

1.5017 

1.5150 

1.2 

0.0822 

0.1251 

0.4481 

0.3922 

0.2044 

0.1014 

0.0403 

— 

1.3937 

1.4416 

1.5 

0.0372 

0.0334 

0.1671 

0.2568 

0.2295 

0.1706 

0.1136 

— 

1.0082 

1.2645 

2.0 , 

0.0105 

0.0068 

0.0394 

0.0837 

0.1089 

0.1165 

0.1042 

— 

0.4700 

1.0036 


(c) Is - 2p eigenstates coupling approximation 


Singlet 

k i 

L = 0 

L = 1 

L - 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2s 

— 

0.9 

1.0 

1.1 

1.2 

1.5 

2.0 

0.0044 

0.0168 

0.0299 

0.0296 

0.0059 

0.0010 

0.1216 

0.0655 

0.0366 

0.0169 

0.0037 

0.0006 

0.1422 

0.3011 

0.3948 

0.3088 

0.0821 

0.0131 

0.0206 

0.0851 

0.1421 

0.0989 

0.0271 

0.0057 

0.0260 

0.0517 

0.0718 

0.0327 

0.0059 

0.0099 

0.0236 

0.0458 

0.0335 

0.0060 

0.0103 

0.0304 

0.0296 

— 

0.2682 

0.4156 

0.5883 

0.5830 

0.3386 

0.1376 


Triplet 

2 S + 2 t 


k i 

L - 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 


0.9 

1.0 

1.1 

1.2 

1.5 

2.0 

0.0002 

0.0016 

0.0037 

0.0055 

0.0059 

0.0028 

0.2066 

0.1078 

0.0540 

0.0249 

0.0027 

0.0002 

0.0005 

0.0020 

0.0060 

0.0098 

0.0133 

0.0073 

0.1651 

0.1599 

0.1446 

0.0868 

0.0311 

0.0363 

0.1002 

0.1336 

0.1231 

0.0559 

0.0187 

0.0365 

0.0695 

0.1100 

0.0700 

0.0217 

0.0341 

0.0849 

0.0772 

— 

0.2073 

0.3315 

0.3820 

0.4220 

0.4267 

0.2445 

0.4755 

0.7471 

0.9703 

1.0050 

0.7652 

0.3820 

0.4867 

0.7508 

0.9836 

1.0529 

1.0215 

0.9156 


(d) Is - 2s - 2p eigenstates coupling approximation 


Singlet | 


L - 0 

L ~ 1 

L = 2 

L= 3 

L = 4 

L= 5 

L = 6 

L = 7 

2 S 

— 

0.9 

1.0 

1.1 

1.2 

1.5 

2.0 

3.0 

4.0 

0.0390 

0.0360 

0.0358 

0.0345 

0.0172 

0.0036 

0.0004 

0.0001 

0.0745 

0.1123 

0.1094 

0.0806 

0.0175 

0.0023 

0.0002 

0.0001 

0.1027 

0.2575 

0.3405 

0.2912 

0.0953 

0.0170 

0.0015 

0.0003 

0.0317 

0.0886 

0.1278 

0.1003 

0.0303 

0.0037 

0.0008 

0.0308 

0.0506 

0.0693 

0.0344 

0.0060 

0.0014 

0.0113 

0.0237 

0.0440 

0.0333 

0.0078 

0.0020 

0.0075 

0.0105 

0.0290 

0.0293 

0.0094 

0.0238 

0.0229 

0.0106 

0.2162 

0.4375 

0.6239 

0.6189 

0.3964 

0.1731 

0.0396 

0.0047 


Triplet 

2 s + 2 t 

q t 

k i 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2 t 

0.9 

1.0 

1.1 

1.2 

1.5 

2.0 

3.0 

4.0 

0.0007 

0.0033 

0.0070 

0.0096 

0.0107 

0.0053 

0.0010 

0.0002 

0.0682 

0.0801 

0.0626 

0.0418 

0.0131 

0.0038 

0.0006 

0.0002 

0.0112 

0.0500 

0.0567 

0.0537 

0.0351 

0.0143 

0.0025 

0.0006 

0.1730 

0.1841 

0.1761 

0.1081 

0.0393 

0.0066 

0.0016 

0.1082 

0.1409 

0.1343 

0.0625 

0.0122 

0.0031 

0.0404 

0.0729 

0.1157 

0.0731 

0.0175 

0.0048 

0.0209 

0.0332 

0.0842 

0.0738 

0.0220 

0.0596 

0.0642 

0.0246 

0.0801 

0.3064 

0.4799 

0.5282 

0.5608 

0.3363 

0.0870 

0.0105 

0.2963 

0.7439 

1.1039 

1.1471 

0.9570 

0.5094 

0.1266 

0.0152 

0.3075 

0.7976 

1.1172 

1.1950 

1.1410 

0.9593 

0.6095 

0.4475 
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1.6 


cases so far considered. For a given total 
angular momentum L , the angular momentum 
of the partial wave which is scattered from 
the 2p state may be L - l , L , and L + 1 . The 
first and the third values correspond to a 
wave function which has the same parity as 
the wave functions in the is and the 2s chan- 
nels; in this case, is even. The 

second value corresponds to a wave function 
with a different parity, and the only process 
that occurs with this parity is the 2p elastic 
scattering; in this case, L - l x - 1 2 is odd. 

We have calculated the 2p-2p cross sections 
for the two cases, and they are listed in 
Table 6. The total cross section is shown 
in Figure 6 (on page 33). Because of the r ~ 2 
potential it is difficult to find the zero energy 
value of this cross section. 

The 2s - 2p transition cross section has application in some plasma and stellar atmosphere 
calculations. The total cross section using the Born approximation is given by Seaton (Reference 
26). In Table 7 (on page 32) we list the partial cross section using the close coupling approxima- 
tion. This table may be found useful in problems in which plasma shielding occurs, where only 
electrons with an impact parameter with a given range can induce the 2s - 2p transition. 

It may be noted that the cross sections for the inverse processes 2s - is , 2p - is , and 2p- 2s 
may be calculated by Equation 35 and the symmetry of the T- matrix. 

In all the tables listed here, k x is the wave number in the is , and k 2 is the wave number in 
the 2s or the 2p channels. The energy, in electron volts, of the incident electron in each channel 
is given by E = 13.6 k 2 , where k could be k x or k 2 . All cross sections are in units of 7ra 0 2 . In the 
different tables, 2 is the sum of the partial cross sections calculated. The total cross section Qj. is 
obtained by adding the contribution of higher partial waves than those calculated using the regular 
Born approximation; this could easily be done with the help of the table of the Born approximation. 



Figure 5 — Is — 2p total excitation cross section. (Is - 2p 
refers to Is - 2p eigenstates coupling approximation; 
Is - 2s - 2p has similar meaning. EX. NEGL. refers to 
exchange neglected case, BORN is the Born approxima- 
tion, and EXP. refers to experiment.) 


CONCLUSION 

The noniterative technique employed here can be applied to a large class of problems containing 
exchange integrals. The method is particularly useful when the exchange potential is comparable to 
the direct potential, in which case the convergence of iteration is slow. 
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Table 6 

2p - 2p Elastic Cross Section. 


(a) L - -i, ~'t 2 odd, Born approximation 


^2 

L = 1 

<N 

II 

■J 

L = 3 

L - 4 

L = 5 

L - 6 

L - 7 

S„ 

0.24 

26.562 

6.1961 

2.1260 

0.92680 

0.44497 

0.22694 

0.08105 

36.56 

0.50 

14.210 

5.2190 

2.1768 

1.0526 

0.55433 

0.31476 

0.18383 

23.71 

0.68 

8.8346 

4.0761 

1.9476 

1.0100 

0.54983 

0.31728 

0.18718 

16.923 

0.83 

6.0065 

3.2062 

1.6993 

0.94122 

0.53204 

0.31331 

0.18680 

12.8854 

1.23 

2.5053 

1.7048 

1.1013 

0.70920 

0.44518 i 

0.28212 

0.17577 

6.9237 

1.80 

0.91169 

0.74615 

0.56949 

0.42475 

0.30158 

0.21122 

0.14245 

3.30733 

2.87 | 

0.23647 i 

0.22647 i 

0.20035 i 

0.17178 

0.13878 

0.10926 

0.08177 

1.16488 


(b) L - -4 2 odd, exchange neglected 2p eigenstates coupling approximations 


k 2 

L = 1 

L = 2 

CO 

ll 

*4 

L - 4 

L = 5 

L = 6 

L = 7 

So 

0.24 

61.12 

8.444 

2.5808 

1.1408 

0.6200 

0.3720 

0.2516 

74.52 

0.50 

15.292 

6.884 

2.5436 

1.1580 

0.6216 

0.3640 

0.2352 

27.10 

0.68 

8.008 

4.940 

2.2548 

1.1160 

0.6160 

0.3640 

0.2352 

17.54 

0.83 

5.108 

3.6364 

1.9292 

1.0380 

0.5968 

0.3592 

0.2340 

12.90 

1.23 

2.0484 

1.7392 

1.1812 

0.7652 

0.4972 

0.3236 

0.2208 

6.776 

1.80 

0.7640 i 

0.7220 

0.5816 

0.4436 

0.3312 

0.2400 

0.1784 

3.261 

2.87 | 

0.2132 | 

0.2180 

0.2004 

0.1776 | 

0.1520 ; 

0.1272 

0.1064 

1.195 


(c) L - odd, 2p eigenstates coupling approximation 


Singlet | 

k 2 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L - 7 

Zos 

0.24 

2.963 

4.161 

0.6725 

0.2861 

0.1552 

0.0930 

0.0629 

8.394 

0.50 

3.735 

3.182 

0.7915 

0.3066 

0.1576 

0.0913 

0.0589 

8.323 

0.68 

2.165 

1.728 

0.6851 

0.3025 

0.1587 

0.0920 

0.0590 

5.190 

0.83 

1.371 

1.107 

0.5562 

0.2799 

0.1547 

0.0913 

0.0589 

3.619 

1.23 

0,5280 

0.4603 

0.3115 

0.1991 

0.1277 

0.0823 

0.0558 

1.7647 

1.80 

0.1928 

0.1831 

0.1475 

0.1124 

0.0837 

0.0605 

0.0449 

0.8249 

2.87 

0.0534 

0.0547 

0.0503 

0.0445 

0.0381 

0.0319 

0.0267 

0.2996 

Triplet 

k 2 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

z or 

0.24 

49.22 

3.265 

1.850 

0.8528 

0.4650 

0.2788 

0.1886 

56.12 

0.50 

7.791 

1.851 

1.504 

0.8193 

0.4598 

0.2722 

0.1764 

12.87 

0.68 

4.720 

2.219 

1.361 

0.7700 

0.4481 

0.2700 

0.1756 

9.964 

0.83 

3.373 

2.099 

1.237 

0.7198 

0.4318 

0.2652 

0.1744 

8.300 

1.23 

1.484 

1.225 

0.8376 

0.5509 

0.3630 

0.2384 

0.1639 

4.863 

1.80 

0.5674 

0.5337 

0.4295 

0.3283 

0.2456 

0.1784 

0.1330 

2.4159 

2.87 

0.1596 

0.1632 

0.1500 

0.1328 

0.1136 

0.0951 

0.0797 

0.8940 


(d) L - ^ even, Born approximation 


k 2 

L - 0 

L = 1 

II 

L=3 

L = 4 

L = 5 

L = 6 

o- 

11 

S E 

Qt 

0.24 

12.488 

230.42 

4.0427 

1.9014 

0.30797 

0.57480 

0.29675 

0.06128 

250.09 

290.69 

0.50 

0.1758 

74.475 

7.0963 

0.92172 

0.48094 

0.45745 

0.30098 

0.05695 

83.965 

109.714 

0.68 

0.07386 

36.478 

6.5573 

1.0156 

0.29902 

0.27895 

0.23113 

0.04364 

44.978 

' 63.742 

0.83 

0.22032 

21.559 

5.4143 

1.1653 

0.26124 

0.15835 

0.15225 

0.02533 

28.956 

43.572 

1.23 

0.25024 

7.3349 

2.9703 

1.1453 

0.37750 

0.11091 

0.04377 

0.03572 

12.269 

20.571 

1.80 

0.13477 

2.3518 

1.3155 

0.74676 

0.39290 

0.18754 

0.07731 

0.02824 

5.2348 

9.6066 

2.87 

0.04395 

0.55809 

0.40443 

0.30267 

0.21866 

0.15188 

0.09753 

0.05832 

1.8355 

3.8272 

3.91 

0.01868 

0.20876 

0.17048 

0.14307 

0.11656 

0.09255 

0.06896 

0.04864 

0.8677 

2.1455 


(e) L 

i 

1 

even, exchange neglected Is - 2s - 2p eigenstates coupling approximation 

k 2 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2 e 

q t 

0.24 

31.96 

91.21 

154.6 

92.86 

55.86 

36.79 

24.37 

— 

487.7 

566.32 

0.50 

9.371 

12.65 

30.88 

20.53 

12.99 

8.451 

5.695 

— 

100.6 

129.8 

0.68 

4.156 

5.953 

13.19 

9.628 

6.563 

4.508 

3.130 

— 

47.13 

66.56 

0.83 

2.542 

4.560 

7.424 

5.249 

3.811 

2.793 

2.032 

— 

28.41 

43.07 

1.23 

1.208 

2.879 

2.906 

1.734 

1.106 

0.8438 

0.6912 

— 

11.37 

19.56 

1.80 

0.5612 

1.385 

1.238 

0.8299 

0.5109 

0.3005 

0.2022 

— 

5.028 

9.382 
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Table 6 (Concluded) 

2p - 2p Elastic Cross Section. 


(f) L --C. -l 2 even. Is - 2p eigenstates coupling approximation 


Singlet 

k 2 

L = 0 

L = 1 

L- 2 

L = 3 

L = 4 

L = 5 

L = 6 

t- 

n 

Z ES 

— 

0.24 

1.964 

5.238 

17.34 

1.260 

0.4896 

0.2460 

— 

— 

26.54 

— 

0.50 

0.5131 

3.159 

1.544 

1.439 

0.4063 

0.1906 

0.1133 

— 

7.365 


0.68 

0.2346 

2.091 

1.039 

0.4944 

0.2539 

0.1520 

0.0984 

— 

4.363 


0.83 

0.1227 

1.252 

0.9930 

0.1532 

0.1191 

0.0990 

0.0766 

— 

2.816 


1.23 

0.1071 

0.6796 

0.6167 

0.2249 

0.0635 

0.0249 

0.0217 

— 

1.7384 


1.80 

0.0571 

0.3366 

0.2926 

0.1818 

0.0954 

0.0426 

0.0172 

— 

1.0233 


Triplet 

k 2 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L= 6 

L - 7 

Zet 

Qt 

0.24 

15.62 

26.72 

32.44 

6.674 

1.252 

0.7359 

— 

— 

83.44 

178.89 

0.50 

4.844 

3.384 

19.97 

12.52 

0.8285 

0.3266 

0.2754 

— 

42.15 

72.80 

0.68 

2.247 

3.886 

10.33 

5.922 

1.035 

0.2270 

0.1764 

— 

23.82 

45.22 

0.83 

1.293 

3.650 

6.287 

3.476 

0.9248 

0.2097 

0.1162 

— 

15.957 

32.45 

1.23 

0.4842 

2.287 

2.424 

1.418 

0.6072 

0.2079 

0.0717 

— 

7.500 

17.28 

1.80 

0.1909 

1.063 

0.9673 

0.6659 

0.4001 

0.2024 

0.0962 

-- 

3.586 

8.943 


(g) L - Z 1 - even, Is - 2s - 2p eigenstates coupling approximation 


Singlet | 

k 2 

L = 0 

L - 1 

L =* 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

^ES 

— 

0.24 

7.852 

13.45 

38.56 

21.41 

15.44 

8.610 

6.432 

— 

111.75 

— 

0.50 

2.470 

5.026 

7.433 

4.900 

3.201 

2.101 

1.417 

— 

26.548 


0.68 

1.344 

2.283 

3.025 

2.683 

1.756 

1.166 

0.7903 

— 

13.047 


0.83 

0.7424 

1.316 

1.579 

1.518 

1.079 

0.7580 

0.5313 

— 

7.524 


1.23 

0.2813 

0.6752 

0.6254 

0.3663 

0.2726 

0.2289 

0.1916 

— 

2.641 


1.80 

0.1357 

0.3339 

0.2862 

0.1816 

0.1065 

0.0659 

0.0471 

0.0383 

1.1952 


2.87 

0.0513 

0.1053 

0.0930 

0.0739 

0.0553 

0.0391 

0.0268 

0.0181 

0.4628 


3.91j 

0.0267 

0.0451 

0.0422 

0.0373 

0.0318 

0.0261 

0.0212 

0.0166 

0.2470 


Triplet 

k 2 

L = 0 

L= 1 

L- 2 

— 

L = 3 

L - 4 

L - 5 

L = 6 

L = 7 

V 

Qt 

0.24 

27.90 

75.79 

63.10 

87.21 

41.31 

26.12 

17.13 

— 

338. §6 

518.92 

0.50 

4.337 

4.540 

20.66 

21.62 

10.38 

6.373 

4.255 

— 

72.17 

122.01 

0.68 

2.823 

4.018 

10.75 

8.986 

5.149 

3.337 

2.302 

— 

37.37 

67.46 

0.83 

2.032 

3.674 

6.400 

4.796 

2.950 

2.013 

1.468 

— 

23.333 

44.54 

1.23 

0.9741 

2.290 

2.465 

1.599 

0.9568 

0.6421 

0.4924 

— 

9.419 

20.11 

1.80 

0.4291 

1.055 

0.9303 

0.6614 

0.4150 

0.2540 

0.1639 

0.1180 

4.027 

9.531 

2.87 

0.1564 

0.3217 

0.2907 

0.2330 

0.1780 

0.1291 

0.0907 

0.0623 

1.4619 

3.945 

3.91 

0.0805 

0.1366 

0.1285 

0.1143 

0.0982 

0.0810 

0.0655 

0.0524 

0.7570 

— 
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Table 7 

2s - 2p Excitation Cross Section. 


(a) Born approximation 


k 7 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

2 

*Qt 

0.245 

210.45 

465.14 

449.77 

343.68 

148.60 

209.74 

158.31 

1985.69 

13560 

0.500 

6.2469 

30.729 

57.561 

64.948 

54.469 

49.590 

40.391 

303.935 

3465.0 

0.678 

0.92003 

6.6619 

16.736 

24.356 

25.702 

24.598 

21.122 

120.096 

1930.9 

0.831 

0.24580 

2.1480 

6.4066 

11.011 

13.641 

14.135 

12.973 

60.560 

1308.4 

1.225 

0.02039 

0.20887 

0.78437 

1.7481 

2.7859 

3.6015 

3.9501 

13.099 

620.51 

1.803 

0.00187 

0.01900 

0.07823 

0.20201 , 

0.38645 

0.59807 

0.78354 

2.0692 

294.95 

2.872 

0.00011 

0.00107 

0.00444 

0.01211 

0.02545 

0.04450 

0.06688 

0.15456 

120.26 

3.905 

0.00002 

0.00017 

0.00067 

0.00182 ! 

0.00389 

0.00699 

0.01096 

0.02452 

66.509 


(b) Exchange neglected Is - 2s - 2p eigenstates coupling approximation 


k. 

o 

n 

■J 

L = 1 

L = 2 

L - 3 

L = 4 

L = 5 

L = 6 

2 

Qt 

0.245 

5.311 

12.59 

23.21 

— 

— 

— 

— 

41.11 

12476 

0.500 

0.8651 

10.55 

1.143 

18.92 

26.25 

— 

— 

57.73 

3308.8 

0.678 

1.150 

5.907 

0.5760 

7.168 

12.28 

14.46 

15.08 

56.62 

1867.4 

0.831 

1.249 

2.859 

0.3831 

3.158 

6.664 

8.591 

9.391 

32.295 

1280.1 

0.225 

0.3847 

0.4266 

0.0881 

0.4560 

1.403 1 

2.328 

3.022 

8.108 

615.52 

1.803 

0.0654 

0.0516 

0.0142 

0.0553 

0.1846 

0.3930 

0.6183 

1.3824 

294.26 


(c) Is - 2s “ 2p eigenstates coupling approximation 
Singlet 


k 2 

t- 1 

ii 

o 

rH 

II 

►4 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

2 S 

— 

0.245 

2.243 

4.424 

3.276 

— 

— 

— 

— 

9.943 

— 

0.500 

0.1241 

1.605 

1.348 

6.360 

7.159 

— 

— 

16.596 


0.678 

0.0362 

1.446 

0.5518 

3.056 

3.693 

3.911 

— 

12.694 


0.831 

0.1866 

0.9881 

0.2615 

1.488 

2.144 

2.433 

2.505 

10.006 


1.225 

0.1048 

0.1516 

0.0384 

0.1639 

0.4584 

0.7031 

0.8529 

2.473 


1.803 

0.0175 

0.0157 

0.0046 

0.0148 

0.0530 

0.1108 

0.1726 

0.3890 


2.872 

0.0014 

0.0010 

0.0004 

0.0010 

0.0032 

0.0075 

0.0133 

0.0278 


3.905 

0.0002 

0.0002 j 

0.0001 | 

0.0002 | 

0.0005 j 

0.0012 

0.0019 | 

0.0043 



Triplet 

2 S + 


k 2 

o 

ii 

h4 

L = 1 

L - 2 

L = 3 

L = 4 

L - 5 

L = 6 

Z T 

0.245 

0.0000 

10.40 

56.74 

— 


— 

— 

67.14 

77.083 

12512 

0.500 

2.322 

7.125 

2.357 

3.838 

17.01 

— 

— 

32.652 

49.248 

3300.3 

0.678 

1.590 

2.363 

0.3333 

1.442 

6.810 

9.868 

— 

22.406 

35.100 

1867.0 

0.831 

0.9885 

1.125 

0.1226 

0.8444 

3.411 

5.518 

6.556 

18.566 

28.572 

1276.4 

1.225 

0.2648 

0.2173 

0.0402 

0.2297 

0.7623 

1.423 

1.987 

4.924 

7.397 

614.8 

1.803 

0.0455 

0.0310 

0.0094 ! 

0.0383 ! 

0.1223 1 

0.2544 | 

0.4070 | 

0.9079 ; 

1.297 

294.18 

2.872 

0.0040 

0.0024 

0.0009 

0.0031 

0.0095 

0.0208 | 

0.0359 1 

0.0766 

0.1044 

120.21 

3.905 

0.0007 

0.0004 

0.0002 

0.0005 , 

0.0015 i 

0.0035 ; 

0.0063 

0.0131 

0.0174 

66.502 
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[14.8451 -ft + ink 2 2 ] 


1 

1_ 7 7 


5 

V” 1 


I— l 

n- 1 


1 

2 In V 


= 1 + 4k V 


(Reference 26). 
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Appendix A 

The Four Differential Equations 


+ k 2. L ( L+ n +J. 

dr 2 1 r 2 r . 


u(kj L, r) 


= 2y 0 (Is Is, r)u(k, L, r) + 2y 0 ( Is 2s , r) u (k 2 L, r) 


r l i i/2 

[_3( 2L '+ i">J Vj (Is 2p, r) u(k 2 L- 1, r) 

r l + 1 i 1/2 

[3(2LTT)J y, (Is 2p, r) u(k 2 L+ 1, r) 


+ H7+T rR io< r )y L (is k, l, r) + a7+T rR 2 o( r )y L ( lsk 2 L - r ) 


2/3 


+ 2/3 


3L 


- 1 1/2 


[_(2L + 1) (2L- l) s 


(r) y L -i (lskj L- 1, r) 


2/S 


3(L + 1 ) ~] 1/2 


(2L+ 1) (2L+ 3) 2 J 
/35(L, 0) (l +kj 2 ) (ls|k 1 L)rR 10 (r) 


rR 2 l( r )yL+i( lsk 2 L+ r ) 


/SS(L, 0) (l +k 2 2 ) (ls|k 2 L)rR 20 (r) 


- /SS(L - 1, 0) (l +k 2 2 ) (ls|k 2 L- l) rR 21 (r) 


(Al) 
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-_df_ 

.dr 2 




u(k 2 L, rj 


2y 0 (Is 2s, r)u(k 1 L, r) +2y 0 (2s 2 s, r ) u(k 2 L, r) 


+ 2 


3(2L+ 1) 


1/2 


yj ( 2 s 2 p, r) u(k 2 L - 1, r) 


L + 1 


|_3( 2L + 1) 

2/5 


1/2 


yj (2s 2p, r) u(k 2 L+ 1 , r) 


2/5 


2 L+T rR io ( r > y L ( 2 s k ! L, r) + 2 lTT rR 2o( r )y L ( 2 sk 2 L, r) 


+ 2/5 


3L 


(2L+ 1) (2L- l ) 2 


rR 


21 ( r ) y L -i (2sk 2 L 1, rj 


2 /5 


3(L + 1) 


1/2 


( 2 L + 1 ) ( 2 L + 3) 2 J 
^S(L, 0) (" 5 ' + k, 2 ) ( 2 s | kj L) rK 10 (r) 

/3S(L, 0)(^-+k 2 2 ) (2s|k 2 L) rR 20 (r) 


rR 2 i (O y L + i (2s k 2 L + 1, r) 


/3S(L - 1 , 0) (- 3 - + k 2 2 ) ( 2 s|k 2 L - 1 ) rR 21 (r) 


(A2) 
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[^2 + k 2 2 ' (L r 2 1)L + r ] U ( k 2 L - !• r ) 


2 [3(2L+1)] y,(ls2p, r)u(k l L, r) + 2 [ 3 ( 2 ^ 1)] y x (2s 2p, r) 
+ 2 j^y 0 (2p2p, r) + y 2 (2p2p, r)J u(k 2 L- 1, r) 


6 VL(L + 1) , 

T " 2L+ 1 y 2 ( 2 P 2 P> r ) u(k 2 L+ 1, r) 


[ 3L 1 1/2 

(2L+ 1) (2L- 1 )2 J rR io( r )y L -i ( 2 P k l L - r ) 

[ 3L “| 1/2 

(2L+ 1) (2L - 1 ) 2 J rR 20 ( r ) y L-i (2pkjL, r) 


6/3 

+ 2L-1 rR 2l( r > 


y L (2pk 2 L- 1, r) L - 1 , 

(2L+1) 2 + 371^-2 ( 2 P k 2 L - 1- r ) 


6/3W.(L + 1) . , 

(2L+ 1)2 Vl ( 2 Pk 2 L+ 1, r) rR 21 ( r ) 


- /8S(L. 1)(-J- + kj 2 ) (2pikj L) rR 10 (r) 

- /3S (L, l)(-J- + k 2 2 ) ( 2p| k 2 L) rR 20 ( r ) 

- /38(L, 2)(i + k 2 2 ) (2p|k 2 L- l)rR 21 (r) 


( k 2 L * r ) 


(A3) 
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+ T u(k 2 L + i, r 


f L + 1 ] 1/2 , f L + 1 1 ,/2 

2 |_3( 2L + 1)J yj (Is 2p. r)u(k,L, r) - 2 [p-qj J y 2 (2s 2p, r) u(k 2 L, r) 


6 r L(L + 1) , , , 

‘ 5 ' 2L + 1 1 yi ( 2 P 2p, r) u (k 2 L - 1 , r) 


(L + 2) y 2 (2p2p, r)" 


2 |_y 0 (2p 2p. r ) + 5 ( 2L + 1 ) J u ( k 2 L+ !. r ) 


, [ 3(L + 1) 1 " 2 

|_(2L + 1) (2L + 3) 2 . 


rR io ( r ) y L +i ( 2 P k i L, r) 


[ 3(L + 1 ) 1 1/2 

(2L+ 1) (2L + 3) 2 rR 2 o( r >y L+1 ( 2 P^L. r ) 


- (2L+ 1)2 rR 2l( r )y L (2 P k 2 L-l. r) 


6/3 [y L (2p k 2 L+ 1, r) (L + 2) y L+2 (2p k 2 L + 1 , r) 

+ 2L + 3 rR 2i ( r ) _ (2L + l) 2 + (2L + 5) 


- /3S(L, 0)1-4 + k 2 2 ) ( 2 P 1 k 2 L + l) r R 2 i ( r ) * 
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Appendix B 

Elements of the Potential Matrix 



Elements of Fjj 


2L+ 1 


"£io f r 

- ' L Jo 


R 10 r ' L+1 dr'-R 10 r 


L+l 


r rR io ; 

Jo^ dr . 
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X. 
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R 24 [ R 20 R 2l] ’ 


3/3 


L(L + 1 ) 
( 2 L + 1 ) 4 


®21 f r /* r ^21 

— R 2 i r' L+1 dr - R 21 r L+1 -J 
_ J o Jo r 


dr ' 


R 34 = *43 . 


Elements of gjj and h|j 


/*!„ r L+ ’ 

811 r 2 L +1 ’ h n 


1 1+k i 2 


~2 S (L. 0) r 


822 


833 


5R r L+1 
“ k 20 r 


2 L + 1 ’ ^22 R 20 

3/3 R 21 r L+1 


T +k 2 2 

2 S ( L - °) r 


( 2 L- 1 ) ( 2 L+ l) 2 


h 33 


21 


3 / 3 (L- 1 )R 21 r L_1 

833 " (2L - 1 ) ( 2 L - 3 ) ’ h 33 _ R 2i 


1 7 + k 2 2 
ir~2 2— *( L - 2 > r 


844 


84$ 


3/3 Rjj r L+1 


4 +k 2 2 


( 2 L + 3 ) (2L+ l) 2 

h 44 R 21 |^ r L 2 

- S(L, 0) rj 

3 / 3 (L + 2 ) R 21 r 1 -^ 

2 _ R 21 


( 2 L + 3 ) ( 2 L+ 5 ) ’ 

h 44 r L+2 



812 2L + 1 R 2o fL * ’ ^12 


- 1 1 + k 2 2 

~^L 2 S < L - °> r 


821 812 [ R 20 ^ R 10 ] * ^21 ^12 [ R 10 ^ R 20 ] 


813 


/ 3/3 


L 

1/2 

p _L 

U - p 

1 l +k 2 2 1 

_( 2L + 1) ( 2 L- l) 2 

K 21 r 

h l 3 K 10 

r L-i 2 5 ( L » 1 ) r 


814 


813 [ R 21 * R lo] ’ ^31 ^13 [ R 10 ^ R 2l] 

L + 1 "I 1 '' 2 


y 5/3 


(2L + 1 ) ( 2 L+ 3 ) 2 


R 2ir L+2 • h l4 = ^Tl 
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&41 g 14 [ R 21 ^ R lo] * ^41 ^14 [ R 


10 * R 2l] * 


g 23 


= ftP 


g 24 


L 

1/2 

T? 


- 1? 

1 ? +k 2 2 

_(2L + 1) (2L- l) 2 

K 21 f ’ 

h 23 

K 20 

TiTi 2 8 (L. 1) r 


g 32 g 23 [ R 21 ^ R 2 o] ’ ^32 ^23 [ R 20 ^ R 2l] 

- f3fi 


L + 1 


( 2 L+ 1 ) ( 2 L + 3 ) 2 


R ai' L+2 » h 24 = 


g 42 

g 24 [ R 21 ^ 

' R 2 o] 

’ ^42 

^24 [ R 20 

* R 2 , 

g 34 

= - 3/3 

1 1 

* 1 ) " 
l ) 4 

1/2 

R r L+1 
*21 r 

• h 34 = 

R 21 

"rL 


g 43 834 


In the F. . matrix the interchange of the functions R 10 , R 20 , and R 21 accompanies the inter- 
change of their arguments too. 


42 



Appendix C 

Removal of the Singularity in the Determinant of the Coefficients 
of the Linear Transformation for L = 0, 1 

Case I: L = 0 

By using the definition of D 4j and F u and Equations 63 in the text, the following relation can 
be derived from Equation 59 (text): 


where the superscript /x is suppressed when there is only one value for u , and 


I 


rR 20^-^7 +k l 2 ) V 1 “Z 3 rR i 


— + k 2 
dr 2 2 


<3r [ a I3 B 24 @ a 23 B mJ • (Cl) 


a l3 = f R 10 R 21 r3dr ° [2 15 X 3-’] ^ * 

/•oo 

a 2 3 ~ ] ^20^21 ° 3 >5 . 

J 0 

Integrating the left-hand side of Equation Cl by parts, and making use of Equations 11 and 63 (text), 
we obtain 



2[B 21 -/3B u ] 


(C2) 


We conclude that 


B 2 i - /3B 12 



- ^ a 23 B I#] 


(C3) 


Equation C 3 connects the right-hand sides of the four equations of Equation 62 specified by 
ij =21, 12, 24, 14. A similar relation should hold among the left-hand sides of these equations. 
This in fact is the case and, by making use of the first of Equations 63, it can be shown directly that 
equations similar to Equation C3 hold among the elements of each column kl^ of the left-hand side 
of Equation 62 specified by ij =21, 12, 24, 14. We conclude that one equation of Equation 62 is 
linearly dependent on others and the determinant of Equation 62 is singular. 
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Case II: L = 1 

Similar to the previous case, the following relation can be derived from Equation 59 in the 

text: 


rI?21 Ur 2 +k ‘ 2 r 2 / ^ r ®l° (dr 2 + ) V 3 


dr 


& a i3 B 11 + & a 23 B 12 3 13 B 33 + ^( 3 13 B 34 


(C4) 


where 


33 


fa 

J o 


4 dr = 30 . 


Integrating the left-hand side of Equation C4 by parts, and making use of Equations 11 and 63, we 
obtain 



„ / d 2 

2 \ ( d 2 ,\ 


r M^ +k *- 

— jv, -/3rR 10 (^— + k/jv 3 

•*0 




dr = - 2[B J1 -/3B l3 ] . 


Combining Equations C4 and C5, we get 

B 31 ^ B 13 ~ 3 j^'' a 13 B 11 

Finally, Equation 59 (text) gives the following relation: 


+ ^ a 23 B 12 3 13 B 33 + ^( a i3 B 34 5 ^ 3 33 B l 


(C5) 


(C6) 


rR 2i( dr 2 + k 2 2 r 2^ V 2 ' /3rR 20^ dr 2 + k 2 2 j V 3 


dr 


3 ^ a 23 B 22 + ^ a l3 B 21 a 23®33 + ^( 3 23 B 34 3’^ 3 33 B 24) 


Integration by parts of the left-hand side gives, as before, 


/*0C 




( d 2 , , 2 \ 

( d 2 ,\ 


21 (dr 2 2 ~ r 2 / 

1 v 2 -/3rR 20 ^-^ + k 2 2 J v 3 

^0 

L_ ' ' 

_J 


(C7) 


dr - 2 [B 32 /3B 23 ] , ( C g) 
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whereupon we get 


®32 ^®23 3” a 23 B 22 + ^ a l3 ®21 a 23®33 + ^( 3 23 B 34 "S’^ a 33®24^ - (C9) 

Similar to the case L = o, Equations C6 and C9 indicate that two equations of Equation 62 are 
linearly dependent on others and the determinant of Equation 62 is singular. 

To remove the singularity in the L = o case, one of the is chosen to be arbitrary; and a 
degenerate equation is removed from Equation 62 (text). Similarly, in the L = l case two of the 
C^j are chosen arbitrary; and two degenerate equations are removed from Equation 62. 
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Appendix D 


Elements of the Matrix of the Sum of the Asymptotic 
Coulomb and Centrifugal Potentials 


U n = L(L + 1 ) r' 2 , = L(L + 1 ) r~ 2 , 


= (L- 

l)Lr' 2 

+ 12(L 

- D(2L + 

l)" 1 r' 3 , 

U 44 = 

(H 

l)(L + 2)r' 2 + 12(L+ 2) (2L+ l)' 1 r’ 3 

= U 21 


• ”i3 


T 256 /2 

L 

1/2 

~ 2 


= 0 

= U 31 = 

[ 2~4r~ 

(2L + 1>_ 


r , 


= U 41 


256/2 

" (L + 1)" 

1/2 

2 

It — 

U 32 

- _ zr 

r l i» /2 _ 2 


243 

_c 2L + n 

r , 

U 23 

0 

_(2L+ 1)J r • 

= U 42 

= 6 

ikLLL 

(2L + 1) 

1/2 

r~ 2 , 

ii 

$ 

»43 = ' 

36 

[L(L + 1)] 

1/2 (2L+ l)' 1 r‘ 3 . 
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